Incidence of nonextensive thermodynamics in 
temporal scaling at Feigenbaum points 



A. Robledo 

Instituto de Fisica, 
Universidad Nacional Autonoma de Mexico, 
Apartado Postal 20-364, Mexico 01000 D.F., Mexico. 



Abstract 

Recently, in Phys. Rev. Lett. 95, 140601 (2005), P. Grassberger addresses the inter- 
esting issue of the appUcabihty of (^-statistics to the renowned Feigenbaum attractor. 
He concludes there is no genuine connection between the dynamics at the critical 
attractor and the generalized statistics and argues against its usefulness and cor- 
rectness. Yet, several points are not in line with our current knowledge, nor are his 
interpretations. We refer here only to the dynamics on the attractor to point out 
that a correct reading of recent developments invalidates his basic claim. 
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1 Introduction 

A sharp disagreement and refutation has been published recently [1] on the 
contention, based on many studies [2] - [23], that the generalized statistical- 
mechanical scheme known as nonextensive statistical mechanics [24] is per- 
tinent to the dynamical properties at the so-called onset of chaos in dissi- 
pative low-dimensional iterated maps. The 'dialogue' between advocates and 
opponents of the appropriateness of this new formalism, here referred to as 
g-statistics, to this old problem has been hampered by the lack of both: i) 
an all-inclusive explanation of how and why the new formalism relates to the 
aforementioned dynamics, and ii) an objective appraisal of the published re- 
sults on the subject as they have developed in time and in depth. Whereas a 
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definite rationalization may be hastened by the ongoing accumulation of in- 
formation from current research, unbalanced judgments are sustained by the 
prevalent belief that, after the historic and extensive work on the transitions 
to chaos carried out a few decades ago, there are no major unknown features 
left to be revealed, nor intrinsic value in new approaches. Here we attempt to 
break this impasse by explaining as plainly as possible the manner in which q- 
statistics manifests at the transitions to chaos, and by indicating how it relates 
to previous approaches, such as the 'thermodynamic formalism' for nonlinear 
dynamics as adapted by Mori and colleagues [25] to the study of this kind of 
attr actor. 

When the control parameter of a nonlinear one- dimensional map is set at the 
threshold between periodic and chaotic motion an unusual and complicated 
dynamics arises, whose features have been explored long ago [26]- [29] and 
described recently in greater detail [16] - [23]. For both the period doubhng 
and the quasiperiodic routes to chaos the transition from periodic to chaotic 
behavior is mediated by the appearance of a multifractal 'critical' attractor 
with geometrical properties already known for a couple of decades [30]-[32]. 
Critical attractors have a vanishing Lyapunov coefficient Ai and a sensitivity 
to initial conditions that does not converge to any single-valued function 
but instead displays a fluctuating pattern that grows as a power law in time 
t [25]. Trajectories within such a critical attractor show self-similar temporal 
structures, they preserve memory of their previous locations and do not have 
the mixing property of truly chaotic trajectories [25]. A special version of the 
thermodynamic formalism for deterministic chaos [31], [33]-[37] was adapted 
long ago [25] to study the dynamics at critical attractors and quantitative 
results were obtained that provided a flrst understanding, particularly about 
the envelope of the fluctuating and the occurrence of a so-called 'g-phase' 
dynamical phase transition [25]. 

Below we describe the scaling properties of the fluctuating sensitivity ^t(xj„) 
at a multifractal critical attractor as well as those of its associated spectra of 
generalized Lyapunov coefficients A(xj„) (where Xin stands for the initial posi- 
tion of trajectories). We illustrate these properties and their connection to the 
g-statistical expressions by considering mainly the specific case of the period 
doubling onset of chaos, i.e. the so-called Feigenbaum attractor. We explain 
that the dynamics at the attractor consists of families of dynamical g-phase 
transitions and discuss the ensuing relationship between the thermodynamic 
and the g-statistical formalisms. 

We refer to various points raised in Ref. [1] as a means of clarifying our cur- 
rent understanding of this problem [38]. We chose this format in part because 
subsequent critical opinions [39] on our developments [16] - [23] appear to be 
based on a cursory and unquestioning belief that the arguments in Ref. [1] 
are indisputable and its views definitive. However, as we demonstrate below. 
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the contents in Ref. [1] suffer shortcomings that seem to stem from overcon- 
fidence and poor reading of the studies cited there, particularly [16] - [23]. 
Another reason for adopting an extended 'comment' layout is that the cus- 
tomary channels for response have been discouraged or barred by the editors 
of the journal and magazine where the criticisms have been published, pre- 
venting an informed, in-depth and balanced debate on this issue. We would 
wish that this propagation of swiftly-taken and disingenuous judgements be 
counterbalanced by standard objective reasoning. 

In Section 2 we sum up the features of g-statistics in relation to the dynamics 
we discuss. In Section 3 we explain the structure of for the fluctuating dy- 
namics at the Feigenbaum attractor. We point out there that this dynamical 
structure is distinguished by a built-in 'aging' or 'waiting time' scaling law. In 
Section 4 we describe the relationship between the thermodynamic approach 
and g-statistics as applied to critical attractors. We indicate that the flxed 
values for the TsaUis' q index correspond to the values that Mori's held vari- 
able q takes at the dynamical phase transitions, i.e. c\trans = Q- In Section 5 we 
demonstrate the condition for the linear growth of the g-entropy Sq, i.e. the 
occurrence of g-phase transitions. We also prove the identity between the rates 
of g-entropy growth and the g-generalized Lyapunov coefficients. For clarity of 
presentation the previous sections focus on the properties that link only the 
two dominant scaling regions of the multifractal attractor, its most crowded 
and most sparse. In Section 6 all the scaling regions of the multifractal are 
considered, it is shown how is obtained in terms of the discontinuities of 
Feigenbaum's trajectory scaling function a, and it is observed that the entire 
dynamics is made of a hierarchical family of pairs of g-phase transitions. In 
Section 7 we refer to the parallelisms in the expression of g-statistics in the 
other two routes to chaos displayed by low-dimensional maps, via quasiperi- 
odicity and via intermittency. In Section 8 we summarize the main results. 



2 g-statistics for critical attractors 



First of all it is indispensable to stipulate what (according to the author) 
is intended by the manifestation of g-statistics in the dynamics at a critical 
attractor. This can be summarized in two related properties. The first concerns 
the (finite time) sensitivity to initial conditions ^t, defined as 



^t{xin) = lim 

Axo—*0 



Axt 



Aa;,v 



(1) 
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where Axj^ is the initial separation of two trajectories and Axj that at time 
t. This quantity is expected to involve expressions of the form 

6(a^m) = exp [\q{Xin) t], (2) 



where q is the entropic index, \q{xin) is the g-generalized Lyapunov coefficient, 
and expg{x) = [1 — (g — l)x]~^^^'^~^^ is the ^'-exponential function. The second 
property relates to 'temporal extensivity' of entropy production [40] at critical 
attractors [2], [20], [22], [23]. That is, linear growth with time t of the entropy 
associated to an ensemble of trajectories. The fitting expression for the rate 
of entropy production Kg{xin) to be used is thought to be given by 

Kq{Xin) t — Sq(t,Xin) — Sq(0, Xin)i (3) 



where 



1 _ q 

Sq^Y^P^^nqpr^^^^, (4) 



is the Tsallis entropy [24] and where Pi{t) is the distribution of the trajectories 
in the ensemble at time t given that they were distributed initially within a 
small interval around Recall that In^y = — 1)/(1 — g) is the inverse 
of expg{y) and notice the explicit dependence of all quantities on the initial 
position Xin- 

Several comments are in order: 

i) It is anticipated that the index q in Eq. (2) takes a well defined value 
determined by the basic attractor properties such as its universal constants. 
For a multifractal critical attractor there may be a distinct discrete family of 
such values for q. 

ii) Because of the memory retention of the trajectories the dependence on Xin 
in Eq. (2) (or Eq. (3)) does not disappear for sufficiently large t. Notice that 
the usual large t limiting condition is not present in the definition of in Eq. 
(2) nor in Eq. (3). 

iii) Since fluctuates according to a deterministic pattern (given by the spe- 
cific route to chaos studied) the variable t in Eq. (2) cannot run through all 
positive integers sequentially but only through specific infinite subsets. How- 
ever, and as seen below, all times t appear in Eq. (2) when Xj„ is varied. 

iv) The standard exponential form for with the ordinary Ai is obtained from 
Eq. (2) when g — > 1. This limit is brought about when the control parameter 
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in the map is shifted from its value at the chaos threshold to a value that 
corresponds to either a periodic or a chaotic attractor. For these cases the 
fluctuations in die out and the Lyapunov spectra Xq{xin) collapses into the 
single number Ai independent of Xin for t large. 

v) Depending on the value of q and the sign of \q{xin), the g-exponential sen- 
sitivity can grow or decrease asymptotically as a power law with t, but it can 
also grow or decrease faster than an ordinary exponential. As we see below, for 
a multifractal critical attractor Eq. (2) is obtained from an exact transforma- 
tion of a pure power law [19], [20], but a faster than exponential g-sensitivity 
is rigorously obtained for the critical attractor at a tangent bifurcation [18]. 

vi) The value of the index q that makes the rate Kg a time independent 
constant is expected to be the same value appearing in Eq. (2). 

vii) The identity K(,{xin) = \q{xin) is supposed to hold. In the limit q ^ 1 
the familiar identity Ki = Ai [32] is recovered (where the rate of entropy 
production Ki is given by Kit — Si{t) — Si{0) and = — Z^iPilnpj). Note 
that in this identity the rate Ki is not the so-called KS entropy [30]- [32] but 
a closely related quantity. See Section 5. 



3 Fluctuating dynamics at the Feigenbaum attractor 

Let us consider now properties specific to the Feigenbaum map g{x), obtained 
from the fixed point equation g{x) = ag{g{x/a)) with g{0) = 1 and g'{0) = 0, 
and where a = —2.5029079... is one of Feigenbaum's universal constants [30]. 
For expediency we shall from now on denote the absolute value \a\ by a. 
Numerically, the properties of g{x) can be conveniently obtained from the 
logistic map 



with jj, = = 1.401155189.... A first important issue that merits amplifica- 
tion and explanation lies behind Grassberger's comment that his Eq. (3) in 
[1], here rewritten as 



holds only for special values of time, e.g. t = 2'^ — 1,A; = 0, 1,... when Xm = 1. 
Not quite, there are many other time sequences all of which satisfy Eq. (6) 
exactly [20]. To see this consider the trajectory Xn with initial position Xin ~ 
0, plotted in Fig. 1 in a way that makes evident the structure of the time 



(5) 



^t{Xin) = Oi^ 



(6) 
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sequences and the power law property that their positions share. It can be 
observed there that the positions Xn{k,i) of the sequences with times n{k, I) = 
{21 + 1)2^, each obtained by running through A; = 0, 1, 2, ... for a fixed value 
of I — 0,1,2,... fall along straight diagonal hnes. The starting elements of 
these sequences are the positions X2i+i- To be precise, the entire attractor 
can be decomposed into this family of position sequences as every natural 
number n appears in one of them. Actually, the complete fluctuating sensitivity 
^t{xin = 1) can be decomposed into a family of identical power laws given by 
Eq. (6) with t — n—{2l+l) [20]. The sequence for Z = belongs to the envelope 
of while those for other values of I are shifted rephcas that describe interior 
values of This feature offers an enormous simpliflcation to the study of the 
dynamics on the attractor. Time rescaling of the form t/tu, with = 2/ + 1 
makes all power law sequences collapse into a single one in a manner analogous 
to the so-called 'aging' property of 'glassy' dynamics where tw is a 'waiting 
time' [41], [42]. 
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Fig. 1. Absolute values of positions in logarithmic scales of the first 1000 iterations 
for the trajectory of the logistic map at the onset of chaos /Xoo with initial condition 
Xin = 0. The numbers correspond to iteration times. The power-law decay shared 
by the time sequences mentioned in the text can be clearly appreciated. 



The family of power laws in Eq. (6) can be rewritten [20], with the use of the 
identity 



j. \ In a/ In 2 



2/ + 1 



(7) 
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with t — (21 + 1)2'^ — 2Z — 1, as the family of g-exponentials, 

et(xi„) = expJA«t], (8) 



all with the same value of q. Above, q = 1 — In 2/ In a and A^'^ = (21 + 
l)"MnQ;/ln2 with t = n — 2/ — 1, when Xi„ = 1, or with t = n = (21 + 1)2^ 
when Xin is any of the positions X21+1 in Fig. 1, in all cases running through 
k — 0, 1, ... while I is fixed to a given value I — 0, 1, ... Alternatively, Eq. (8) 
can be written as the single ^'-exponential 

et(a;,„) = expjAf-^]. (9) 



So, contrary to the statement in [1] - that his Eq. (3) is not a scaling law 
because it only holds for special values oit - Eq. (6), and its equivalent forms 
Eqs. (8) or (9), in fact correspond to a more complex form of a scaling law, one 
that fits the fluctuating dynamics at the attractor, and referred here as a 'two- 
time' scaling law. Below wc make clear that the positions of these sequences 
for large k belong to one scaling region within the multifractal, its most sparse 
region. 

A persuasive reason for considering a g-exponential as a further way to write 
a power law is the presence of a time scale factor: the generalized Lyapunov 
coefficient A^'^ that appears multiplying t in This useful quantity (hidden 
in the power law) can be immediately 'read' from the anomalous just as 
the ordinary Lyapunov coefficient Ai is read from the exponential of chaotic 
dynamics. It is important to stress that it is in the scaling limit (large n 
or k) that becomes a'', and it is this which transforms (exactly) into the 
g-exponential Evidently, reference and use of this property involves more 
than a mere choice between a power law and a ^-exponential representation of 
the sensitivity. The beauty of a scale-free power law is there to enjoy. As for 
generality, all of these properties apply with only minor changes to the period- 
doubling transition to chaos in unimodal maps with extremum of order z > 1 
[22]. The derivation of Eq. (6) for z > 1 is given in [21] so there is no need to 
reproduce it in [1]. 

So, what is the advantage in studying the nonmixing trajectories at the Fcigcn- 
baum attractor in terms of the above-explained power-law network structure 
for instead of the time and position averages described in [1]? Clearly, the 
dynamical organization within the attractor is hard to resolve from a simple 
time evolution: starting from an arbitrary position on the attractor and 
recorded at every t. What is observed are strong fluctuations with a scram- 
bled structure that persist in time. Conversely, unsystematic averages over 
Xin and/or t would rub out details of the multiscale properties. However, if 
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specific initial positions with known location within the multifractal are cho- 
sen, and subsequent positions are observed only at pre-selectcd times, e.g. 
t = 2^ — 1, when the trajectories visit another region of choice, a well-defined 
g-exponential sensitivity appears, with q and the Lyapunov spectrum X{xin) 
fixed by the attractor universal constants. 



4 Thermodynamic approach and g-statistics 



We now address the connection between the thermodynamic approach as em- 
ployed by Mori and colleagues [25] to study dynamical phase transitions at 
critical attractors and the g-statistical properties of the Feigcnbaum attrac- 
tor. In particular: i) We establish the equivalence of the generalized Lyapunov 
coefficients present in the two methods, ii) We show the link between the oc- 
currence of dynamical phase transitions and the existence of unique values for 
the entropic index q. iii) We demonstrate the linear growth of the Tsallis en- 
tropy at the dynamical phase transitions and the identity between the entropy 
growth rates with their corresponding g-Lyapunov coefficients. 



4-1 Generalized Lyapunov coefficients and q- deformation 



We make use of the scaling features of just described to examine Grass- 
berger's advice in [1] that Eq. (4) in [1] gives the natural generalization of the 
Lyapunov coefficient Ai. Indeed, a closely related prescription. 



(10) 



was initially adopted in Refs. [27] , [29] . It differs however from Eq. (4) in [1] in 

that no time average is taken. As a function of consecutive values of t, \t{xin) 
does not converge to any constant but does undergo recurrence. This recur- 
rence suggests inspection of Xt{xin) for values of t along the time sequences 
described in the previous section. A straightforward calculation shows that 
Xt{xin — 1) as defined in Eq. (10) is exactly given by 



Xt{Xin = 1) = jlnq 



dXir, 



In a 



{21 + l)ln2 



A(0 



(11) 



where t, as before, runs through the sequences t = (21+1)2'^— 21— 1, k = 0,1, 
with / = 0, 1, ... fixed. So, the earlier definition for the generalized Lyapunov 
coefficient is equivalent to that given for the same quantity by the g-statistics. 
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The meaning of the index q is given above. It is the degree of 'g-dcformation' 
of the ordinary logarithm that makes Aj finite for large t. It is not difficult to 
corroborate that this result is valid for general on the attractor. 



^.2 Dynamical partition function 



Eqs. (10) and (11) suggest a broader connection between the thermodynamic 
formalism and the g-statistical features found for the dynamics at the Feigen- 
baum attractor. Indeed, it has been shown in [22], [23] that the physical origin 
of the existence of a distinct value (or values) for the entropic index q is linked 
to the occurrence of dynamical phase transitions, of Mori's g-phase type [25], 
which connect qualitatively different regions of the attractor. 

The thermodynamic formalism is built on the dynamical partition function 
Z(n,q) = j d\ P{X,t) W{X,ty-'^, (12) 



where the weight W{X, t) was assumed [25] to have the form W{X, t) = exp(At) 
for chaotic attractors and W{X, t) — t^ for the attractor at the onset of chaos, 
with A a generalized Lyapunov coefficient (whose dependence on t and Xin is 
not written explicitly). The variable q plays the role of a 'thermodynamic field' 
like the external magnetic field in a thermal magnet or, alternatively, 1 — q 
can be thought analogous to the inverse temperature. We note that W{X,t) 
can be written in both cases as 



W{X,t) 



(13) 



where f{x) is the iterated map under consideration, that is. 



Xfix, 



Bit) In 



(14) 



with B{t) = t^^ for chaotic attractors and B{t) = (Int)^^ for the fiuctuating 
sensitivity at the critical attractor [25]. As mentioned, for chaotic attractors 
the fluctuations of Xt{xin) die out as t ^ oo and this quantity becomes the 
ordinary Lyapunov coefficient Ai independent of the initial position Xi^. For 
critical attractors Xt{xin) maintains its dependence on t and Xin but as al- 
ready seen it is a well-defined constant provided t takes values along each time 
sequence, or if the two-time scaling property (Eq. (9)) is summoned. 
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For critical attractors the density distribution for the values oi X, t ^ 1, 
P{X,t), is written in the form [25] 



(15) 



where ip{X) is a concave spectrum of the fluctuations of A with minimum 
■^(0) = and is obtained as the Legendre transform of the 'free energy' func- 
tion 0(q), defined as 0(q) = — limt^oo In q)/ Int. The generalized Lya- 
punov coefficient A(q) is given by A(q) = d(f){q)/dq [25]. The functions 0(q) 
and il^{X) are the dynamic counterparts of the Renyi dimensions -D(q) and the 
spectrum f{a) that characterize the geometric structure of the attractor [31], 
[32]. For non hyperbolic attractors like those at the onset of chaos the func- 
tions 0(q) and 4'{X) that are used to describe the dynamics are independent 



4-3 q-phase transitions and Tsallis q index 

As with thermal 1st order phase transitions, a g-phase transition is indicated 
by a section of linear slope rric = 1 — g in the spectrum (free energy) ip{X) 
and consequently a discontinuity at q = g in the Lyapunov function (order 
parameter) A(q). For the Feigenbaum attractor a single g-phase transition 
was numerically determined [29] and found to occur approximately at a value 
around rric — —(1 — q) — —0.7. It was pointed out in Ref. [27], [29] that this 
value would actually be nic — —(1 — q) — — In 2/ In a = —0.7555... 

From the knowledge we have gained on g-gcneralized Lyapunov coefficients 
for this attractor, e.g. Eq. (11), we can determine the free energy ■0('^)- But 
to do this we need to extend our results a bit further. We note that Xt{xin) 
in Eq. (11), corresponds to trajectories that originate in the most crowded 
region of the attractor, e.g. Xin = 1, and terminate at its most sparse region, 
Xn={2i+i)2i' — 0, large k. Naturally, these trajectories when observed at times 
n — {21 + 1)2'^ grow apart from each other as k increases {Xt{xin) > 0). 
When the inverse situation is considered, e.g. Xj„ = 0, with final positions 
Xn={2i+i)2k+i - 1, one obtains [22] 



of D(q) and f{a). 



in 



) 



expQ[Ayt], 



(16) 



with Q — 



2 - 



g = l-|-ln2/lnQ; and 



Xt{x, 



0) 



21na 




(17) 



in 



(2Z + l)ln2 
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where t = (2/ + 1)2^ -21, k = 0, 1, / = 0, 1, ... Above, of course, At(0) < 
since the trajectories when observed at times (2/ + 1)2^ + 1 come progres- 
sively close to each other. Notice that the change in sign of and the 
relation Q — 2 — q match the property of inversion of the g-exponential, 
expg{x) — 1/ exp2_^(— x). Also, trivially, when the initial and final positions 
of trajectories belong to the same region within the attractor, e.g. xi — 1 and 
Xn={2i+i)2k+i ^ 1, one has = 1 with = 0. 

From the results for A^'^ and Xq the Lyapunov function A(q) and spectrum 
■0(A) are constructed [22] 



A(q) = 



A(o), -oo < q < g, 
0, q<q<Q, 

Xq\ Q < q < oo. 



(18) 



and 



1 - Q)X, Ag) < A < 0, 
(l-g)A, 0< A< Ajo). 



(19) 



with A(°) = In a/ In 2 ~ 1.323 and X^q' = -2Xf^ ~ -2.646. The constant 
slopes of ip{X) represent the g-phase transitions associated to trajectories link- 
ing two regions of the attractor, x ~ 1 and x ~ 0, and their values 1—q and q—1 
correspond to the Tsallis index q obtained from C,t- The slope g — 1 ~ —0.7555 
coincides with that initially detected in Refs. [28], [27]. 

The significance of the above exercise is that it reveals the physical reason for 
the existence of a well-defined value of the entropic index q in critical attractor 
dynamics. This is the occurrence of a g-phase transition. The value of the field 
q at which the transition takes place is precisely q (and the same of course for 
the inverse process, C{trans = Q)- 



5 Temporal extensivity of the g-entropy 



As we have seen, along each time sequence n{k,l) = {21 + 1)2^, I fixed, the 
sensitivity with Xin = 1 is given by the g-exponential in Eq. (8), and con- 
sequently we can express Xt{xin) in terms of the g-deformed logarithm in Eq. 
(11). Therefore for times of the form t — n — {21 we have 

P{\t)W{\t) = 5{X - A«) exp^(A(')i), (20) 
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and 



Zit,q) = W{X^^\tY-^ 



(1 - <l)^^t 



(l-q)/(l-q) 



(21) 



Notice that in Eq. (21) q is a running variable while g = 1 — In 2/ In o; is fixed. 

Our next step is to consider the uniform probability distribution, for fixed A^'^ 

and t, given by pi{\^^\t) — W ^ , i — 1, ...,W , where W is the integer nearest 
to (a large) W . A trajectory starting at = 1 visits a new site Xn belonging to 
the time sequence n{k, I), I fixed, every time the variable k increases by one unit 
and never repeats one. Each time this event occurs phase space is progressively 
covered with an interval of length Axk^i — Xi^2i+i)2'' ~ ^(21+1)2''+^ > 0. The 
difference of the logarithms of the times between any two such consecutive 
events is the constant In 2, n large, whereas the logarithm of the corresponding 
phase space distance covered Axkj is the constant In a. Wc define Wk to be the 
total phase space distance covered by these events up to time n = {21 + 1)2^. 
In logarithmic scales this is a linear growth process in space and time from 
which the constant probability Wj, of the uniform distribution pj is defined. 
In ordinary phase space x and time t we obtain instead Wf^ — exp{—X'^H). 
We then have that 

w 

Z{t,q) = EbiW]' = 1 + (1 - c\)S,{t), (22) 

i=l 



where S'q = luq W is the Tsallis entropy for pj. 

The usefulness of the g-statistical approach is now evident when we recall 
from our discussion above that the dynamics on the critical attractor displays 
as its main feature a g-phase transition and that at this transition the field 
variable q takes the specific value qtrans — Q- Comparison of Eqs. (21) and 
(22) indicates that at each such transition both Z and Sq grow linearly with 
time along the sequences n(k,l). In addition, A^^^ can be determined from 
Sq/t = . Thus, with the knowledge we have gained, a convenient proce- 
dure for determining the relevant quantities for the fiuctuating dynamics at 
a critical multifractal attractor could take advantage of the above properties. 
Curiously, the same features have been observed for Sq with q — Q associated 
to trajectories (with linear ^t) in a conservative two-dimensional map with 
vanishing ordinary Lyapunov coefficients [43]. 

Further, the identity Sq/t — A^'^ (also derived in [22]) between the rate of 
Q'-entropy change and the generalized A is not the identity 



r^(5i(i) -5i(o)) = rMn 



dxi^ 



(23) 
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in [1] (the zero identity for t — > oo) but refers to Xt{xin) as above. Yes, it 
considers an instantaneous entropy rate but it is comparable in the sense of 
[44] to the g-generahzed KS entropy studied in [12]. The identity Sq/t = A^'-* 
holds for t — > oo as the interval length (around Xin) vanishes. It fluctuates, 
but as explained, we look for the detailed dependence on both Xj„ and t. The 
invariant density is not smooth and is built up by placing an initial condition 
at each point Xin on the attractor. In contrast to the chaotic case there is 
not one identity but many, and the claim in [1] that averages are needed for 
applications of Pesin's identity seems ineffectual for nonmixing trajectories. 
Our results may not be insignificant as these coincide (when I — 0) with those 
in [12] where the q-KS entropy was considered. On the contrary, the entropies 
in [1] from symbolic dynamics do not sense (i.e. do not depend on) the 
universal a and/or the nonlinearity z. 



6 Hiercirchical family of gr-phase transitions 

With reference to the 'rich zoo' of values for the index q alluded in [1] , there 
is a well-defined family [22] of these within the attractor that deserves to be 
explained. This family is determined by the discontinuities of Feigenbaum's 
trajectory scaling function a (which measures the convergence of positions 
in period 2^ orbits as k ^ oo) [30] and are all expressed in terms of the 
universal constants of unimodal maps (for general nonlinearity z > 1). Eqs. 
(8) and (16) and the values of q and 2 — q given there are obtained from the 
largest discontinuity in a, that is itself related to the most crowded and most 
sparse regions in the attractor. The other discontinuities lead to expressions 
for similar to Eqs. (8) and (16). There is a corresponding family of pairs of 
Mori's g-phase transitions, each associated to orbits with common starting and 
finishing positions at specific locations of the attractor. As before, the special 
values for q in are equal to those of the field variable q in Mori's formalism 
at which the transitions occur [22] . Since the amplitudes of the discontinuities 
of a diminish rapidly, there is a hierarchical structure and consideration only 
of the dominant discontinuity provides a sound description of the dynamics. 
When all discontinuities are considered it becomes apparent that the dynamics 
on the critical attractor is constituted in its entirety by the infinite but discrete 
set of gr- phase transitions. We reproduce below some basic expressions [22]. 

The trajectory scahng function a is obtained as the A; — > oo hmit of 

^.(J) = (24) 

where the numbers d^j are the so-called diameters that measure the bifurca- 
tion forks that form the period-doubling cascade sequence [30] . The diameters 
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arc determined from the positions of the 'superstable' periodic orbits of lengths 
2'^, i.e. the 2''-cycles that contain the point a; = at /I^ < fJ'oo [30]. This func- 
tion has finite (jump) discontinuities at all rationals of the form yj = j/2^~^^, 
and we denote the values of a before and after them (omitting the subindex 



k) as cr{y~) — l/a~ and (y{y^) — '^/oij, respectively. The numbers af are 
universal constants. For the largest discontinuity yj — one has — o? and 
OL^ — —a. 

The strategy employed [22] in determining from a is to chose the initial and 
the final separation of the trajectories to be the diameters Axj^ = d^^j and 
Axj = 4j+t, t^2^ - I, respectively. So, Ct(xj„), Xj„ = Xinij), is obtained as 



it{xin) = lim 



d. 



k,j+t 



d, 



k,j 



(25) 



Notice that in the /c — > oo limit Axin 0, t ^ oo and the 2^-supercycle 
becomes the onset of chaos (the 2°°-supercycle). Then, for each discontinuity 
of a at yj, Ctixin) can be written as [22] 



(26) 



For the inverse process, starting at Aa;j„ = dkj+t = —dk,j-i and ending at 
Axt = dkj = —dkj-i+t, with t = 2*^ + 1 one obtains [22] 



(^kiVj ) 



, ^ = 2^= + !, A;> 1 



(27) 



Therefore 



a 



+ 



6,: 



f28) 



and, similarly to Eq. (6), the sensitivities in Eq. (28) can be re-expressed as 
g-exponentials with g-indexes 



q=\ 



In 2 



In 



, g=2-g= 1+ 



In 2 



In 



(29) 



respectively. All other consequences described in the previous sections follow 
in a straightforward manner. 
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It is interesting to note that discrete infinite families of values for the index 
q arise in different contexts, such as in the proposed g-generalization of the 
ordinary and Levy-Gnedenko central limit theorems [45]. 

As a check on the leading role of the most crowded and sparse regions of 
the attractor in determining its dynamics we consider a two-scale Cantor set 
approximation of the multifractal attractor of fractal dimension df. Following 
a standard procedure [31] for a trajectory xq, Xi, a;2fc+i_i, k large, on the 
attractor we cover it with a set of intervals of lengths ll = \xi — a;j_,_2fe|, 
i — 0, ...,2'^ — 1. The smallest length scale Z^^^ ~ a~'^^ is observed in the 
neighborhood of a; = 1 while the largest ~ occurs near x = Q. 
Trajectories starting in the vicinity of x = 1 in the multifractal visit the 
vicinity of a; = at times of the form t = 2'^ — 1. In the two-scale set the 
equivalent trajectories of duration t = 2'^ — 1 starting at positions assigned to 
the scale l^^^ visit positions that correspond to the scale /^2x ^ number of times 
equal to dft. For these times \dxi/dxQ\ — l^l^/l^in otherwise \dxi/dxo\ — 1. 
Use of this in the time-averaged expansion rate [1] 



1 



t=2'^-l 



dxo 



k > 1 (30) 



leads to 



m 



Plnt^ dfln^,k^l (31) 



or 



p = dflna/ln2. (32) 

As df ~ 0.5388 [30] we obtain /3 ~ 0.7131 which falls (given the number of 
digits considered in the value of df) within 0.033% of the numerical value 
reported for /3 in [1]. Eq. (32) is indeed a very simple relation linking the 
expansion rate constant (3 and the universal constant a. 



7 g-statistics and the other routes to chaos 



With reference to the generality of our results for the Feigenbaum attractor, 
and besides all unimodal maps of arbitrary nonlinearity z > 1, there is a very 
similar picture obtained for another important multifractal critical attractor. 
This corresponds to the quasiperiodic route to chaos, recently studied [23] in 
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the framework of the famihar golden-mean (gm) onset of chaos in the critical 
circle map [30], [32]. The dynamics on the attractor (a fat fractal [46]) is more 
involved than that for the Feigenbaum case but there is a strong parallelism 
that can be glimpsed through the following equivalences when only the most 
crowded and most sparse attractor regions are considered [23]: 

i) The time sequences along which the sensitivity to initial conditions exhibits 
the power law scaling ^t{xin = 1) = a'^ are n = (/ — m)Fk + 'mFk-2, where 
Fk is the Fibonacci number of order k, each sequence is obtained by running 
through k — 1,2, 3, ... for fixed values of Z = 1, 2, 3, ... and m = 0, 1, 2, I — 1. 
Notice that the sequences now depend on one additional index. As before there 
is a time shift t = n — t^, here with t.^, = / — m + mWg^, linking the two scales 
t and n and with the property that rescaling of the form makes all power 
law sequences collapse into a single one. The number a is now the absolute 
value of the universal constant agm — —1.288575 obtained as the scaling factor 
that satisfies the fixed-point map equation g{x) — agmg{o(gm9{x / 0(lm) ■ 

ii) The role of the period doubling time scale factor 2^ is now taken (asymp- 
totically) by where w"^ = (\/5 — l)/2 is the golden mean. 

iii) The sensitivity is given by g-exponentials as in Eqs. (8) and (16) with 
q = 1 + lnwg„i/2\na , A^^'"*^ = 2 In «/((/ — m) -|- mWg^) Inwgm, Q = 2 — q and 

Ag''"^ = -2 A J'-). 

iv) A pair of ^-phase transitions take place at q = g and at q — Q — 2 — q, 
that correspond to switching starting and finishing orbital positions. And, as 
before, linear growth (or reduction) of the Tsallis entropy occurs when 
q = 5 (or q = Q) with rates i^f - AJ'-'") (or K^q'""^ = Ag'""^). 

More generally, use of the trajectory scaling function a for the golden mean 
quasiperiodic attractor reveals an infinite family of pairs of g'-phase transitions, 

with a hierarchical structure, each member with properties like in i-iv above 
[23] . As before the dominant behavior arises from the most crowded and most 
sparse regions of the multifractal. 

Furthermore, it is significant to call attention to the fact that the fixed-point 
map solution of g(x) = ag{g{x/a)) for the tangent bifurcations of unimodal 
maps of general nonlinearity z > 1, the third route to chaos, is rigorously 
given by a g-exponential map. See [17], [18]. This feature leads immediately 
to a ^'-exponential with = 3/2 for all z > 1 [18]. The tangent bifurcations 
display weak insensitivity to initial conditions, i.e. power-law convergence of 
orbits when at the left-hand side {x < x^) of the point of tangency x^. However 
at the right-hand side {x > Xc) of the bifurcation there is a 'super-strong' sensi- 
tivity to initial conditions, i.e. a sensitivity that grows faster than exponential 
[18]. The two different behaviors can be couched as a g-phase transition with 
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indexes q and 2 — q for the two sides of the tangency point. Also a two time 
or aging scaling property like that in Eq. (9) holds for this critical attractor 
[47] . The sensitivity is dependent on the initial position Xin or, equivalently, 
on its waiting time t^', the closer Xj„ is to the point of tangency the longer 
tw but the sensitivity of all trajectories fall on the same g-exponential curve 
when plotted against t/t^. It is worthwhile to underscore that in this case the 
g-exponential is not an alternative way to express a power law but the exact 
function that describes the faster than exponential increase of 



8 Conclusions 

Thus, we have explained in detail the dynamics at the Feigenbaum and other 
critical attr actors. In all cases the fluctuating sensitivity to initial conditions 
has the form of infinitely many interlaced g-exponentials that fold into a sin- 
gle one with use of the two-time scaling t/ty, property. More precisely, there 
is a hierarchy of such families of interlaced g-exponentials; an intricate (and 
previously unknown) state of affairs that befits the rich scaling features of a 
multifractal attractor. Therefore, the comment in the abstract of Ref. [1] that 
"the... behavior at the Feigenbaum point based on non-extensive thermody- 
namics... can be easily deduced from well-known properties of the Feigenbaum 
attractor" is somewhat an overstatement. 

Summing up, the incidence of g-statistics in the temporal scaling at the Feigen- 
baum attractor is verified; the g-exponential sensitivity is explained; the origin 

of the index g is made clear; and the equality between g-generalized Lyapunov 
coefficients and g-entropy growth rates is demonstrated. A similar corrobora- 
tion applies to the other two known routes to chaos, those via quasiperiodicity 
and via intermittency. 

Acknowledgments. Partially supported by CONACyT and DGAPA-UNAM, 
Mexican agencies. 



References 

[1] P. Grassberger, Phys. Rev. Lett. 95, 140601 (2005). 

[2] C. Tsallis, A.R. Plastino and W.-M. Zheng, Chaos, Solitons and Fractals 8, 885 
(1997). 

[3] U.M.S. Costa., M.L. Lyra, A.R. Plastino and C. TsaUis, Phys. Rev. E 56, 245 
(1997). 



17 



[4] M.L. Lyra and C. Tsallis, Phys. Rev. Lett. 80, 53 (1998). 

[5] F.A.B.F. de Moura, U. Tirnakli and M.L. Lyra, Phys. Rev. E 62, 6361 (2000). 

[6] V. Latora, M. Baranger, A. Rapisarda and C. Tsallis, Phys. Lett. A 273, 97 
(2000). 

[7] U. Tirnakh, G.F.J. Ahahos and C. Tsallis, Phys. Lett. A 289, 51 (2001). 

[8] E.P. Borges, C. Tsahis, G.F.J. Ahahos and P.M.C. de Oliveira, Phys. Rev. Lett. 
89, 254103 (2002). 

[9] G.F.J. Ahahos and C. Tsallis, Phys. Rev. Lett. 93, 020601 (2004). 

[10] J. Yang and P. Grigolini, Phys. Lett. A 263, 323 (1999). 

[11] M. Buiatti, P. Grigolini and L. Palateha, Physica A 268, 214 (1999). 

[12] S. Montangero, L. Franzoni and P. Grigolini, Phys. Lett. A 285, 81 (2001). 

[13] M. Coraddu, F. Melloni, G. Mezzorani and R. Tonelh, Physica A340, 234 (2004). 

[14] M. Coraddu, M. Lissia and R. Tonelh, Physica A 365, 252 (2006). 

[15] R. Tonelh and M. Coraddu, Eur. Phys. J. B 50, 355 (2006). 

[16] A. Robledo, Europhys. News 36, 214 (2005). 

[17] A. Robledo, Physica D 193, 153 (2004). 

[18] F. Baldovin and A. Robledo, Europhys. Lett. 60, 518 (2002). 

[19] F. Baldovin and A. Robledo, Phys. Rev. E 66, 045104(R) (2002). 

[20] F. Baldovin and A. Robledo, Phys. Rev. E 69, 045202(R) (2004). 

[21] E. Mayoral and A. Robledo, Physica A 340, 219 (2004). 

[22] E. Mayoral and A. Robledo, Phys. Rev. E 72, 026209 (2005). 

[23] H. Hernandez-Saldaha and A. Robledo, Physica A, in press, and 
cond-ma t/0507624 

[24] C. Tsallis, J. Stat. Phys. 52, 479 (1988). 

[25] H. Mori, H. Hata, T. Horita and T. Kobayashi, Prog. Theor. Phys. Suppl. 99, 
1 (1989). 

[26] P. Grassberger and M. Scheunert, J. Stat. Phys. 26, 697 (1981). 
[27] G. Anania and A. Politi, Europhys. Lett. 7, 119 (1988). 

[28] T. Horita, H. Hata, H. Mori and K. Tomita, Prog. Theor. Phys. 81, 1073 (1989). 

[29] H. Hata, T. Horita and H. Mori, Prog. Theor. Phys. 82, 897 (1989). 

[30] H.G. Schuster, Deterministic Chaos. An Introduction, 2nd Revised Edition 
(VCH Publishers, Weinheim, 1988). 



18 



[31] C. Beck and F. Schlogl, Thermodynamics of Chaotic Systems (Cambridge 
University Press, UK, 1993). 

[32] R.C. Hilborn, Chaos and nonlinear dynamics, 2nd Revised Edition (Oxford 
University Press, New York, 2000). 

[33] T.C. Halsey, M.H. Jensen, L.P. Kadanoff, I. Procaccia, and B.I. Shraiman, Phys. 
Rev. A 33, 1141 (1986) and references therein. 

[34] J.-P. Eckmann and I. Procaccia, Phys. Rev. A 34, 659 (1986). 

[35] M. Sano and Y. Sawada, Prog. Theor. Phys. 76, 945 (1986). 

[36] H. Fujisaka and M. Inoue, Prog. Theor. Phys. 77, 1334 (1987). 

[37] G. Paladin and A. Vulpiani, Phys. Rep. 156, 147 (1987). 

[38] This paper is a considerably 

expanded version of A. Robledo, cond-mat/0510293 , For other comments on 
Ref. [1] see C. Tsallis, cond-mat/0511213, 

[39] R. Balian and M. Nauenberg, Letter to the Editors, Europhys. News 37, 9 
(2006). Briefly, the studies described in Ref. [16] refute the letter's claim that 
" the value of q is never selected through theoretical arguments" . The verbatim 
quotation "...obey the features of the g-entropy" in point 2 is wrongly copied. 
The fact that the word "chimera" appears nowhere in Ref. [1], in contrast with 
the statement contained in this letter, is symbolic of the shallow and fabricated 
quality of the criticism. 

[40] S. Abe and Y. Nakada, Temporal extensivity of Tsallis entropy and the bound 
on entropy production rate, |cond-mat / 0603550] 

[41] A. Robledo, Phys. Lett. A 328, 467 (2004). 

[42] F. Baldovin and A. Robledo, Phys. Rev. E 72, 066213 (2005). 

[43] G. Casati, C. Tsalhs and F. Baldovin, Europhys. Lett. 72, 355 (2005). 

[44] V. Latora, M. Baranger, Phys. Rev. Lett. 82, 520 (1999). 

[45] S. Umarov, C. Tsallis and S. Steinberg, cond- mat/060 3593 ; S. Um arov, C. 
TsaUis, M. Gell-Mann and S. Steinberg, cond-ma t /0606038| [cond^Imit; /0606040l 

[46] J.D. Farmer, Phys. Rev. Lett. 55, 351 (1985). 

[47] A. Robledo, Molec. Phys. 103, 3025 (2005). 



19 



